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Introduction 

O , After about three decades of extensive efforts in both nuclear experiments and theoretical 

calculations, the equation of state (EOS) of isospin symmetric matter is well understood 

through experiments of collective flow and subthreshold kaon production. Nowadays, the 
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pq ■ nuclear EOS of asymmetric nuclear matter has attracted a lot of attention. The EOS of 

asymmetric nuclear matter can be described approximately by 

O . E{p, S) = Eo{p, 6 = 0) + E,y^,{p)6' (1) 

where 6 = ^2_£e jg isospin asymmetry, Eo(p, 6) is the energy of pure symmetric nuclear 
matter, and Esyjn(p) is the symmetry energy, where Esym(po) = 32 MeV is the symmetry 

^ ! 

^ ! energy at normal nuclear matter density. The symmetry energy is E(p,l) - Eo(p,0), i.e., 

00 . 

;0 : the dffiereuce of the energy per .mcleon between pure neutron matter and symmetr.e 

r^ , nuclear matter. The symmetry energy is important not only to the nuclear physics com- 

munity as it sheds light on the structure of radioactive nuclei, reaction dynamics induced 
by rare isotopes, but also to astrophysicists since it acts as a probe for understanding 
the evolution of massive stars and the supernova explosions. The existing and upcoming 
radioactive ion beam (RIB) facilities lead a way in understanding nuclear symmetry en- 
ergy. Experimentally, symmetry energy is not a directly measurable quantity and has to 
be extracted from observables which are related to symmetry energy. Therefore, a crucial 
task is to find such observables which can shed light on symmetry energy. A large num- 
ber of studies on the symmetry energy of nuclear matter have been done during the past 
decade [1]. In this paper, we aim to see the sensitivity of collective transverse in-plane 
flow to symmetry energy and also to see the effect of different density dependencies of 
symmetry energy on the same. The various forms of symmetry energy used in present 
study are: Egym oc Fi{u), Egym oc F2{u), and Egym oc F^^u), where u = — , Fi(u) oc m. 
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F2{u) oc u^'^, F^{u) oc u"^, and F4 represents calculations without symmetry energy. The 
study is carried out with IQMD model [2]. 

Results and discussion 

We simulate several thousands events for the neutron-rich system of ^^Ca+"^^Ca at energies 
of 100, 400, and 800 MeV/nucleon at impact parameter of b/bmax=0.2-0.4. 

In Fig. 1 we display the time evolution of < p'^" > for different symmetry energies used 
in this paper at 100 MeV/nucleon for particles lying in the different bins. To understand 
the sensitivity of transverse momentum to the symmetry energy as well as its density 
dependence in the Fermi energy region, we calculate the transverse flow of particles having 
— < 1 (bin 1) and particles having — > 1 (bin 2). Dotted, solid and dashed lines 
represent all particles, binl and bin2. Panels (a), (b), and (c) are for Egym oc PjP^'"^, and 
p^, respectively. Panel (d) is for calculations without symmetry energy. The total < p'^^^ > 
is negative during the initial stages and continues decreasing till 30 fm/c which indicates 
dominance of attractive interaction. In panels (a) and (b), it becomes positive whereas 
in panels (c) and (d) it remains negative during the course of the reaction. If we look at 
< p'^^^ > of particles lying in bin 1 for Fi(-u) [Fig. 1(a)] and F2{u) [Fig. 1(b)] in the time 
interval to about 20-25 fm/c, we see that it remains positive. It increases with time upto 
15 fm/c and reaches a peak value. This is because in the spectator region (where high 
rapidity particles lie) the repulsive symmetry energy will accelerate the particles away 
from the overlap zone in the transverse direction. After 15 fm/c, < p'^" > (of particles 
in bin 1) begins to decrease. This is because these particles will now be attracted toward 
the central dense zone. From 10 to 20 fm/c particles in bin 2 continue increasing in the 
midrapidity region. In the case of Fi(-u) and F2{u), particles which enter the central dense 
zone (bin 2) have already a high positive value of < p'^" > (i.e., going away from the dense 
zone). So, the attractive mean field has to decelerate the particles first, make them stop, 
and then accelerate the particles back toward the overlap zone. At about 20-25 fm/c 
particles from bin 1 have zero < p'^^^ > [see shaded area in Fig. 1(a) and 1(b)]. Up to 30 
fm/c, particles feel the attractive mean field potential after which the high density phase 
is over; i.e., in case of Fi(-u) and F2{u) between and 30 fm/c, particles from bin 1 are 
accelerated toward the overlap zone only for a short time interval of about 5 fm/c, whereas 
for the case of F^{u) [Fig. 1(c)] and F4 [Fig. 1(d)] between 0-30 fm/c, particles from bin 
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Figure 1: The time evolution of < p'^^^ > for different forms of symmetry energy for 
different bins at b/bmax=0.2-0.4. 



1 are accelerated towards the overlap zone for a longer time interval of about 20 fm/c 
between 10-30 fm/c. Moreover, the < p^"' > of particles lying in the bin 1 [for F^lu) and 
F4] follows a similar trend. This is because for p/po < 1 the strength of symmetry energy 
F3(m) will be small and so there will be less effect of symmetry energy on the particles, 
which is evident from Fig. 1(c) where one sees that the < p'^^^ > remains about zero 
during the initial stages between zero to about 10 fm/c. The < p'^}^ > due to particles in 
bin 2 (dashed line) decreases in a very similar manner for all the four different symmetry 
energies between and 10 fm/c. Between 10 and 25 fm/c, < p'^^^ > for F3(m) and F4 
decreases more sharply as compared to in case of Fi(m) and F2{u). This is because in 
this time interval particles from bin 1 enter into bin 2. As discussed earlier, < p'^" > of 
particles entering bin 2 from bin 1 in the case of Fi(m) and F2(m) will be less negative due 
to stronger effect of symmetry energy as compared to the case of F3(m) and F4 [3]. 
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